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ABSTRACT 

We carry out the first time-dependent numerical MagnetoHydroDynamic modeling of an extrasolar 
planetary system to study the interaction of the stellar magnetic field and wind with the planetary 
magnetosphere and outflow. We base our model on the parameters of the HD 189733 system, which 
harbors a close-in giant planet. Our simulation reveals a highly structured stellar corona characterized 
by sectors with different plasma properties. The star-planet interaction varies in magnitude and 
complexity, depending on the planetary phase, planetary magnetic field strength, and the relative 
orientation of the stellar and planetary fields. It also reveals a long, comet-like tail which is a result 
of the wrapping of the planetary magnetospheric tail by its fast orbital motion. A reconnection 
event occurs at a specific orbital phase, causing mass loss from the planetary magnetosphere that can 
generate a hot spot on the stellar surface. The simulation also shows that the system has sufficient 
energy to produce hot-spots observed in Ca II lines in giant planet hosting stars. However, the short 
duration of the reconnection event suggests that such SPI cannot be observed persistently. 
Subject headings: stars: magnetic field - stars: coronae - planet-star interactions 



1. INTRODUCTION 

The newly discovered exoplanets are not only newly 
discovered worlds literally, but are also new worlds in 
terms of the physical system they introduce us to. Many 
exoplan ets have now been observed since their first dis- 
covery (Mayor fc Queloz 1995 Schneider] 1 1995 Mayor 



et al. 20l)3p . In particular, many Jupiter-like giant plan- 



ets have been observed at a distance of less than 0.1 AU 
from their parent star, and some even within 10i?*, 
where they ar e essentially located inside the exte nded 



stellar corona (Schneider 1995 Mayor et al. 2003). It 



is not unreasonable to suppose that these gi ant plan 



ets have a substantial internal magnet ic field ( Sanchez 
Lavega|[2004| Durand-Manterola 2009) and a magneto 
sphere, which could encompass a significant fraction of 
the extended coronal volume. Indeed, the Jovian mag- 
netosphere is the largest entity in our solar system after 
the solar one ( |Bagenal et al.|2004 |. 

Perhaps the most fundamental difference between the 
close-in planet scenario and magnetospheres in our own 
solar system is that close-in planets and their magneto- 
spheres can be located within the stellar Alfven radius. 
In analogy to the hydrodynamic sonic point (where the 
flow speed equals the sound speed) , the Alfven radius (or 
the Alfvenic point) is the distance at which the accelerat- 
ing stellar wind equals the Alfven speed, u^ = B / ^/Anp. 
Here B is the magnetic field strength and p is the plasma 
mass density. Beyond the Alfven point, the flow is super- 
Alfvenic and magnetic information (or energy) cannot 
propagate back towards the star and affect the corona. 
Close-in planets could then have a substantial effect on 
the structure of the stellar corona, and the electromag- 
netic interaction between the planetary magnetic field 
and the corona might generate some observational signa- 
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tures. 

In recent years, some signatures of star-planet inter- 
action (SPI) have been observed. In a series of papers, 
|Shkolnik et"aT| ( |2003[ |2005a|bl |2008[ ) have described the 
observed modulations in the Ca 11 K emission line, which 
is an indicator of chromospheric activity. They found 
a line intensity increase corresponding to chromospheric 
"hot spots" with a period that is correlated with the 
planetary orbital motion on HD 189733, HD 179949, r 
Boo, and v And. In some cases the location of these 
hot spots was aligned with the star-planet radial vector, 
while in other cases, the location of the intensity peak 
was shifted from this vector by 70-170 degrees. Some of 



these observati ons also revealed an on/off nature (Shkol- 
nik et al.||2008D. 



At higher energies, Kashyap et al. (|2008|) performed a 



statistical survey of systems with close-m giant planets 
which revealed that the X-ray flux of these systems is 
about 30-400% higher than the typical fluxes from sim- 
ilar star s with planets located furt her out in the stellar 
system. |Poppenhaeger et al. ( 2010[ ) did not find a similar 
effect for a sample of nearby stars, though they did un- 
cover weak evidence for a combined effect of planetary 
mass and orbital distance. 



(2008J) observed 
ound that emis- 



Saar et al. 

the HD 179949 system in X-rays and 
sion associated with a hot spot likely engendered by SPI 
contributed «30% to the emission at a mean plasma tem- 
perature of « 1 keV. A similar trend in the UV band has 



been recently presented by (Shkolnik 2010) 



Stellar coronae are dominated by the stellar magnetic 
field, which al so dictates the structure of the stellar wind 
(Parker 1958). The presence of a close-in planet with 
a significant internal magnetic field could perturb this 
stellar field and might affect the large-scale structure of 
the stellar corona. The planet can affect the corona in 
different ways. First, a purely hydrodynamic affect is 
the planet acting as an obstacle, deflecting the flow of 
plasma (i.e., the stellar wind) around it. Second, an 



2 



electrostatic effect is the large-scale potential field topol- 
ogy being modified by the superposition of the stellar 
and planetary magnetic fields. Third, a MagnetoHydro- 
Dynamic (MHD) effect is expected due to the orbital 
motion of the 'external' field (the planetary field) inside 
the background conducting plasma of the stellar corona 
and the stellar wind. In any case, the process involves a 
modification of the coronal magnetized medium via en- 
ergy transfer between the planetary magnetic field and 
the corona. 

Few theoretical SPI models have been developed to 
explain the observed hot spot s, as well as the overall 
increase in observe d X-ray flux ( Cuntz et al. 120001 |Cran- 
mer & Saar||2007|). In particular, |Lanza| ([2008[ |2009|) 



In particular 

proposed a mechanism to explain the observed phase 
shift in the location of the spots. This work also demon- 
strated that energy transfer from the magnetic field to 
the plasma via magnetic reconncction between the coro- 
nal and planetary field (as the planet is moving through 
the corona) can provide the observed overall energy in- 
crease (about 10 21 W). This amount of energy requires 
rather strong stellar and planetary magnetic fields or al- 
ternatively, this energy can be provided assuming the 
planet triggers magnetic reconnection in a stellar corona 
that is already energized due to accumulation of mag- 
netic helicity provided by photospheric motions. 

Performing numerical simulations of SPI is challeng- 
ing. First, one needs to provide a dynamic model for 
the ambient stellar corona and stellar wind. Second, it is 
necessary to introduce the planet, which is an additional 
boundary condition in the model. Third, for capturing 
the dynamics due to the planetary motion, one needs to 
consi der the p l anet a s a time- dependent bo undar y condi 



tion. 



et al 
in the 



Ip et al | ( |2004[ ), |Preusse et al.| ( |2006 
(2006fnave studied the structure o: 



and Preusse 



Alfven wings 

close-in planet around HD 179949 and made the 
analogy of the magnetic interaction between Ip and the 
Jovian magnetospher e. Lipatov et al. (20051 and Jo- 
hansson et al. (2009) have studied the structure of the 
planetary magnetosphere using a hybrid code driven by 
an approximation to the stellar wind as a boundary con- 
dition. In both cases, the simulations included only the 
planetary magnetosp here 

(hereafter Pa 
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In a recent paper, Cohen et al. 
per I) have presented two simplified IViriJj simulations of 
SPI. In the first case ("Case A"), both the stellar and 
planetary magnetic fields were represented by magnetic 
dipoles, and the planet was fixed in the inertial frame 
of reference. A relative motion due to stellar rotation of 
the coronal plasma with respect to the planetary mag- 
netic field was obtained. In the second case ("Case B"), 
the planetary magnetic field was dipolar again, while the 
stellar magnetic field was obtained using high-resolution 
solar magnetic field data. This provided a realistic stellar 
field, in contrast to the idealized dipole field used in Case 
A. In Case B, the planet was fixed the whole time and 
the simulation was calculated in the frame of reference 
that is rotating with the star in a tidally-loc ked manner 



(like the r Boo system, Butler et al. 1997). The syn 



thetic X-ray light curves of Case A revealed that there is 
a drop in the X-ray flux when the planet is behind the 
star. In Case B, the intensity drop was shifted by ps 60° 
from the star-planet plane. The simulation also showed 



an increase in the total X-ray flux compare to a reference 
case without a planet. 

In paper I, we concluded that the increase in X-ray flux 
in systems with close-in planets is due to the planet pre- 
venting the corona from expanding, so that field lines 
that would be open remain closed. As a result, the 
plasma cannot escape and the overall coronal density is 
higher than in the case without a planet, or with a planet 
located further from the star. Since the X-ray flux is es- 
sentially a line of sight integration of the square of the 
(electron) density, the observed X-ray flux is higher in 
such systems. In addition, in paper I we proposed that 
the observed hot-spots are associated with the magnetic 
connectivity between the star and the planet, and that 
the shift in the location of the hot-spots from the star- 
planet vector is probably due to the complexity of the 
stellar magnetic field 



In a follow-up study to Paper I, Cohen et al.| (2010c 
hereafter Paper II) investigated the influence of chang- 
ing the semi-major axis of the co-rotating planet on the 
stellar mass and angular momentum loss. The disrup- 
tion of the corona and wind occurs when the separation 
is sufficiently small that the planetary and stellar Alfven 
surfaces start to interact. We found that the spin-down 
of stars harboring close-in planets is reduced compared 
to that of stars with distant planets or no planets at all. 

While the results of the simulations presented in Pa- 
pers I and II are significant, both cases were station- 
ary and did not capture all of the dynamical effects that 
are predicted in the theoretical descriptions of SPI (see 
above). In particular, magnetic reconnection between 
the stellar field and the coronal field was not well cap- 
tured, even in Paper I Case A where some relative mo- 
tion between the magnetic bodies was introduced. In this 
paper, we present a dynamical MHD simulation, where 
we include time- dependent circular orbital motion of the 
planet around the star. The model is based on the ob- 
served properties of the HD 189733 system. However, we 
emphasize that the main goal of this study is to inves- 
tigate the interaction of a close-in giant planet with the 
stellar corona of its host star as a fundamental plasma 
physics problem, rather than to model HD 189733 specif- 
ically. 

The structure of the paper is as follows. We present 
the numerical model in Section [2] The results are pre- 
sented in Section [3j and our main findings are discussed 
in Section 21 The results are summarized in Section [U 

2. NUMERICAL SIMULATION 

We perform the numerical simulation using the BATS- 
R-US global MHD model, orig inally developed for the 
solar corona and the solar wind dPowell et al.|1999" 



Toth 



|et al.|[2005|JCohen et al.||2007[ Solar Corona module of 
the Space Weather Numerical Framework). This model 
is driven by synoptic maps of the photospheric radial 
magnetic field (magnetograms), which are routinely ob- 
served on the Sun. Similar observations are available for 
some stellar syst ems using the Zeeman-Doppler-Imaging 
(ZDI) techn ique (jD onati fc Collier Cameron 1997 Donati 
|et al.|1999} see { 2T ). Constrained by this boundary con 



lition tor the magnetic field, as well as the stellar surface 
plasma density, po, surface temperature, To, mass, M*, 
radius, i?*, and rotation period, P*, the model solves the 
set of conservation laws for the mass, momentum, mag- 
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netic induction, and energy (the MHD equations). The 
stellar parameters of HD 189733 used in this simulation 
are listed in Table [T] The end result of the simulation is 
a self-consistent solution for the solar/stellar corona and 
solar/stellar wind, where the volumetric heating needed 
for the wind acceleration is specified according to the 
magnetic field distribution and an empirical relation be- 
tween the solar win d terminal speed and the magnetic 



flux t ube expansion (Wang & Sheeley 1990 Arge & Pizzo 
20001. The implementation of the solar model for the 



stellar c oron al case can be f ound in paper I,|Cohen et al.| 



(2010b), and Cohen et al 



(2010a I. We refer the reader 



to these references tor a detailed description of the nu- 
merical approach taken in the simulation described here. 

2.1. THE STELLAR MAGNETIC FIELD 

In paper I, we approximated the stellar magnetic field 
as a d ipole. However, recent observations by |Fares et al.| 
( 2010 ) provide us with a better description of the stellar 
field in the form of stellar magnetograms of HD 189733. 
These measurements were taken during June and August 
2006, June 2007, and July 2008. In this work, our main 
goal is to study the dynamic interaction between the 
planetary magnetosphere and a complex, realistic stellar 
magnetic field as the former passes through different sec- 
tors of the corona exhibiting different magnetic field and 
stellar wind properties. Therefore, we have a dopted mag- 
netogr ams based on the data presented by |Fares et al.| 
(2010). These low-resolution maps do not contain any 



tine structure, so they are sufficient for our goal We base 
our simulations on the synoptic m ap for 2006, since this 
is the epoch when observations by Shkolnik et al. ( 2008 1 
were made. Figure [T] shows the the map used in our sim- 
ulation; this is directly comp arable to th a t illu strated 
in the top-left of Figure 5 of |Fares et al. (2010). Due 
to the inclination of the star, the map is incomplete in 
its southern hemisphere. However, the missing data oc- 
cupies minor parts of the total area, so we choose to 
superimpose the map with a weak background 5 G (po- 
lar) dipole field, similar to that of the Sun, oriented with 
positive at the North pole, to complete the map. 

Figure [2] shows the steady-state MHD solution for the 
adopted magnetogram without the planet. It shows the 
three-dimensional structure of the stellar coronal field, 
as well as the equatorial plane colored with contours of 
number density, n, and the B z component of the field. 
Phase angles of 0, 90, 180, and 270, which refer to the ref- 
erence positions of the planetary orbit, are also displayed. 
In the remainder of this paper, we use this orientation 
to describe a particular phase angle. Since the planetary 
magnetic field is a perfect dipole directed towards the 
North stellar pole, the map of B z shown in Figure [2] is 
used as the stellar reference when describing the inclina- 
tion between the stellar and the planetary field. 

To act as a source of reference in helping to understand 
the role of the complex surface magnetic topology in the 
simulations including the planet, we also computed a 
simulation for a simplified dipolar stellar magnetic field 
configuration, with an equatorial field strength of 5 G. 



2.2. MODELING THE PLANET 

In our simulation, the planet is represented by an addi- 
tional set of boundary conditions for a second body (the 



first body being the star) . In order to include the orbital 
motion in the simulation, the coordinates of this second 
body are updated according to the time-dependent plan- 
etary orbital motion. Here we assume a circular orbit in 
the equatorial plane, but in principle any inclined orbit 
can be specified. At each time step, the grid structure 
is updated to determine which grid cells are boundary 
cells and which cells are real (non-ghost) cells. The fre- 
quency of the updating of coordinates is constrained by 
the numerical stability condition of the simulation (the 
relation between the grid size, Ax, the time step, At, 
and the Alfven speed). In addition, the grid resolution 
around the second body should be high enough so that 
the body is well-resolved. For this reason, we choose a 
planetary radius of R p = 2Rj, with Rj being the radius 
of Jupiter, instead of the observed value of X.lRj. The 
grid size around the second body is Air = 0.025-R* with 
18 grid cells across the body. Using the flexible grid ca- 
pabilities of the model, we also used this high-resolution 
gridding around the stellar surface (where Aa: = 0.04i?*), 
and along the planetary trajectory. Figure [3] shows the 
grid structure on the equatorial plane with the two bod- 
ies marked in red. 

The strength of the planetary field is unknown. On 
the one hand, based on stellar evolution models the in- 
ternal magnetic fields of hot Jupiters are expected to be 
lower than that of Jupiter due to tidal locking and a con - 
sequently longer rotation period ( |Sanchez-Lavega|2004[ ) . 
On the other hand, a strong field may be required to 
shield these gas giants from erosion by the strong stellar 
wind. Here we study two set of parameters for the planet. 
A "Strong" case, where we use a dipole equatorial field 
of 2 G, and a "Weak" case, with a dipole equatorial field 
of 0.1 G. While we use a fixed boundary condition for 
the planeta ry temperature, T p = 10 4 K (Murray-Clay 
et al.|[2009 ), the boundary condition for the surtace den- 
sity enables us to control the planetary thermally driven 
outflow. This outflow is obtained in the simulation due 
to the fact that the second body serves as a permanent 
mass source in an environment with lower density. Here 
we choose a boundary value of n p = 10 9 cm" 3 for the 
"Strong" case and a boundary value of n p = 10 7 cm~ 3 
for the "Weak" case. In Section H] we discuss the impor- 
tance of this free parameter in the solution. 

3. RESULTS 

A salient characteristic of the simulations carried out 
here is the dependence of the results on orbital phase. 
We emphasize two main aspects here. Firstly, we find 
that when the planet is not co-rotating, as was assumed 
in Paper I, it can no longer provide a steady means of 
preventing the wind-driven plasma from expanding out- 
ward. While the planet still acts to inhibit angular mo- 
ment loss, a consequence of the interaction is that this 
varies through the orbital cycle. Secondly, magnetic re- 
connection occurs at certain phases in the combined stel- 
lar and planetary magnetic fields that allows planetary 
gas to escape the form its magnetospheric confinement. 

In Paper I, we considered the case of a planet that 
was tidally-locked — fixed in location relative the stellar 
surface and co-rotating with the stellar magnetic field. 
This led to the stellar wind being permanently obstructed 
in the direction of the planet and the plasma here being 
bottled up. When the planet moves relative to the stellar 
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field, the wind in the direction of the planet does not 
remain bottled up, but can slip past once the planet has 
moved on. Consequently, the amount of material that 
is confined by the presence of the planet is reduced, and 
the structure of the ambient corona and stellar wind is 
quickly recovered. 

In Paper II, we found that the stellar mass and angu- 
lar momentum loss rates are reduced due to the plan- 
etary effect on the stellar corona. Once the planetary 
and stellar Alfven surfaces start to interact, the coro- 
nal and wind structure is disrupted by the obstacle — the 
planet — so that the mass flux carried by the stellar wind 
decreases. This leads to a long-term feedback from the 
planet on the star in the form of a reduction in the stel- 
lar angular momentum loss rate to the wind, so that the 
spin-down of stars harboring close-in planets is not as 
high as it wo uld be without the plan et, consistent with 
observations ( |Pont|[2009 Lanza 2010). This study again 
considered planets co-rotating with the stellar surface. 
In the simulations here, that restraint has been removed. 

3.1. ORBITAL PEA SE-DEPENDENT MA SS L OSS 

The gas number density from the simulation at four 
different orbital phases is illustrated in Figure [4] for the 
"Strong" case. The salient feature of the density dis- 
tribution is the presence of a comet-like tail of plasma 
trailing the planet. This feature has been observed by 
Vidal-Madjar et al. ( 1994 1 a nd has been predicted from 
hydrodynamic simulations by Schneiter et al. ( 2007 ). We 
will discuss this feature in detail in § |3.3| below. In Fig- 
ure [5] we show the total outward mass flux carried by the 
stellar wind as a function of planetary orbit phase angle 
that underlies the behavior seen in Figure [4] We choose 
to calculate the mass flux through three spheres in or- 
der to illuminate the dependence of the mass loss as a 
function of the distance from the star. The first sphere is 
of radius r — , which crosses the region between the 
star and the planet, the second is for r — 10i?*, which 
encompasses the planet but crosses the planetary magne- 
tosphere, and the third is for r = 23i?*, which occupies 
most of the simulation domain. At r = 6R±, most of 
the mass flux is carried out by the stellar wind, with a 
minor contribution from the planet. The mass flux at 
r = 10R* is dominated by the planetary contribution (in 
addition to stellar wind mass source) , and this is the rea- 
son for its relatively high values. However, most of this 
plasma is actually trapped inside the magnetosphere so 
that it does not contribute to the total mass flux of the 
system — the apparent outward flow from the planet here 
is balanced by a similar flow inward. The total outward 
mass flux at r = 23i?* is instead again dominated by the 
stellar wind, except for phase a ngle s of - 120° , where a 
magnetic reconnection event (§ 3.2) allows plasma to es- 
cape from the planetary magnetosphere and to be added 
to the total radial mass flux. Most of this plasma falls 
onto the star, and this can be seen in Figure [4] as a bridge 
of higher density joining the star and planet. However, 
our mass flux calculations indicate that there is also a 
significant planetary contribution from this reconnection 
to the total mass flux of the stellar wind. 

3.2. MA GNE TIC RECONNECTION AND 
PLANETARY PLASMA ESCAPE 



As the planet moves into a region where the stellar 
field has a different structure, reconnection can occur 
between planetary and stellar fields if they have oppo- 
site polarity. In the "Strong" case, this occurs around 
phase angles of 60 - 120°, when the planet gets into a re- 
gion of negative B z stellar field associated with a helmet 
streamer. This is seen in Figure|2j that illustrates the so- 
lution for the stellar field alone and does not include the 
planet. The planetary plasma is generally confined and 
prevented from expanding freely by the magnetic field 
of the planet. The reconnection event releases some of 
this bottled-up plasma and it subsequently falls onto the 
star. This is seen as a high density region between the 
planet and the star in Figure [4j This high density region 
is absent at phase angles 150 - 300°, prior to the recon- 
nection and after the released plasma has dissipated. In 
the "Weak" case, shown in Figure [6j the planetary mag- 
netosphere suffers from stronger variations than in the 
"Strong" case as seen by the change in the size of its 
tail. The highest mass escape from the magnetosphere 
occurs around 130°, not 60° (for the "Strong" case). This 
demonstrates how sensitive the "observable" feature is 
to the planetary parameters. In this particular simula- 
tion, a weaker planetary field with the same stellar field 
resulted in a smaller deviation of the reconnection site 
from the star-planet axis. This is because a weaker plan- 
etary field "feels" the opposite stellar field later than the 
"Strong" case. 

The planetary outflow and the resulting "hot spot" 
where it falls on the stellar surface for the "Strong" case 
can be seen in Figure [7j where a snapshot is presented for 
the time a magnetic reconnection event has been iden- 
tified. A local region of infalling mass is seen clearly in 
the mass flux plot at the stellar surface. For reference, 
we show a similar snapshot from the "idealized" simu- 
lation computed for the dipole stellar magnetic field of 
5 G and a planetary boundary number density value of 
n p = 10 10 cm~ 3 . In this case, there is a continual open- 
ing and closing of field lines as the planet revolves about 
the star such that there is always an open field channel 
joining the star and planet. The planetary gas is free to 
expand along this open field region and results in a con- 
tinuous infall of matter onto a spot on the stellar surface 
that lies at a constant phase offset from the planetary or- 
bital phase. In the lower panels of Figure [7] this flow can 
be seen impacting the star approximately 180° ahead of 
the planet. 

The reconnection event for the realistic stellar magne- 
togram simulation is shown in detail in Figure [HJ The 
four panels on the left show selected magnetic field lines 
and color contours of the number density over the equa- 
torial plane. When the planet is far from the recon- 
nection region (top-left panel, phase angle of 270°), the 
field topology is of the ambient stellar field, with some 
closed field lines of the helmet streamer and some open 
field lines dragged by the stellar wind. At a phase an- 
gle of 60° (top-right), the field lines that were originally 
open are now connecting the star and the planet. In 
addition, some of the closed field lines of the stellar hel- 
met streamer are now open as a result of the reconnec- 
tion in front of the planet. At a phase angle of 100° 
(bottom-left), reconnection occurs behind the planet, 
closing down some field lines. Finally, at a phase angle 
of 180° (bottom-right), the field topology reverts back to 
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its original, ambient stellar field. The two right panels 
show the reconnection point in front of the planet (phase 
angle of 60°), releasing planetary material onto the star, 
and the reconnection point behind the planet (phase an- 
gle of 100°), causing a disconnection of a plasmoid from 
the planetary tail. 

3.3. PLANETARY MAGNETOSPHERE 
CONFIGURATION 

The shape of planetary magnetospheres in our solar 
system is largely defined by solar wind conditions and, 
in particular, by the orientation of the Interplanetary 
Magnetic Field (IMF) . In the case of close-in planets like 
the system simulated here, the magnetospheric shape is 
defined by the plasma conditions in the particular coro- 
nal sector interacting with the planet. The isolation of 
the planet from the harsh coronal environment depends 
on the magnetopause location, which is defined as the 
surface of pressure balance between the stellar wind dy- 
namic pressure and the planetary magnetic pressure. If 
a planet is located far from the star, the magnetopause is 
associated with a shock in front of it as the stellar wind 
is slowed down to subsonic speeds. In the case of a close- 
in planet, the magnetopause can be located in regions 
where the stellar wind never reaches supersonic speed 
and is still accelerating, so the structure of the surface of 
pressure balance can be highly complicated. For display 
purposes, we choose instead to show the surface where 
the radial velocity in the frame of reference centered on 
the planet vanishes in the direction of the star. Alterna- 
tively, this is the location where the stellar wind in the 
direction of the planet vanishes. Figure [9] shows how this 
location changes in four selected frames of the simulation 
(before, during, and after the reconnection event). It can 
be seen that this boundary moves away from the planet 
following reconnection due to the plasma escaping from 
the planetary magnetosphere. This line never reaches too 
close to the planet, where the magnetic pressure domi- 
nates. Due to the dipole-like magnetic field lines close 
to the planetary surface, we conclude that the dynamic 
pressure of the planetary outflow at this region is not sig- 
nificant. The boundary value for the planetary density 
affects mainly the filling time of the planetary magneto- 
sphere after a reconnection event that allows magneto- 
spheric plasma to escape. 

Another important magnetospheric feature obtained 
from our simulation is the long magnetospheric comet 
-like tail that stretches round almost 90° behind the 
planet. This tail is important since it increases the 
coronal volume affected by the planet. In our static, 
tidily-locked simulations, as well as in the interaction 
of solar wind and magnetospheres in the solar system, 
the magnetospheric tail is aligned with the direction 
of the stellar wind. Here, however, it is perpendicular 
to the wind direction. Therefore, as long as magnetic 
reconnection does not occur, it serves as an elongated 
obstacle standing in the way of the stellar wind, and can 
have a larger impact on the wind structure and mass 
loss rate. The convection electric field, E = — u sw x B 
(u sw is the stellar wind velocity) in this case should be 
stronger than when the magnetospheric tail is aligned. 
The magnetospheres of close-in planets should, then, 
introduce different configurations and dynamics in terms 
of polar cup potentials, particle drifts, and magne- 



tospheric/stellar wind drag. The interaction of such 
magntospheres with Coronal Mass Ejections (CMEs) 
should also be quite different from the interaction of 
CMEs with magnetospheres in the solar system, since 
the CME here will hit the magnetosphere from its side 
and not from its nose. These processes are potentially of 
considerable interest, but are beyond the scope of this 
paper. 

4. DISCUSSION 

Our simulation reveals a highly dynamic interaction 
between the stellar and planetary field as the planet is 
orbiting the star. This is due to the complex stellar field, 
as well as the stellar wind topology. The simulation 
demonstrates that, even when we use a low-resolution 
map for the surface magnetic field, which does not in- 
clude fine structure and small-scale active regions, the 
stellar corona contains sectors with quite different plasma 
properties. The interaction between the planetary and 
coronal field changes from one sector to the next, while 
some of the notable features are due to the planetary 
transition between the different sectors. This is in con- 
trast to the behavior of the "idealized" stellar dipolar 
field simulation, in which the interaction between fields 
is, not surprisingly, more continuous in nature, with an 
open field channel always joining the two bodies. 

The dynamical, short-lived nature of the different fea- 
tures of the interaction in the realistic field case suggests 
that searching for observable SPI signatures could be 
problematic, and that "static" theoretical models for the 
origin of these signatures might not be valid. Based on 
the simulation, the interaction between the coronal field 
and the planetary magnetosphere will not in general be 
persistent, at least for systems with short planetary or- 
bital periods. These are the cases where SPI is expected 
to be the strongest. In addition, the stellar magnetic 
field topology determines the location where changes in 
the interaction take place, but this magnetic topology is 
expected to exhibit secular change within time-scales of 
months or even less. 

The idealized scenario in which the stellar and plane- 
tary fields are described as perfect, opposite dipoles can 
explain persistent particle acceleration as a result of mag- 



netic reconnection between the two fields (Bp et al. 
ICuntz et al.| |2000| ICranmer & Saar||2007T 



2004 
2008 



asma 



|2009[ ). A tidaliy locked system can produce a pi 
density enhancement in closed coronal loops that would 
be open without the planet (paper I). However, based 
on our simulation, none of these proposed scenarios can 
be persistent for the case of a close-in planet that is 
not tidally-lock ed with the stellar rotation. Moreover, 
recent papers (|Zarka| [2007| 



Vidotto et all 120101 have 



suggested that reconnection events between stellar and 
planetary fields can generate radio emissions that can 
be observed and constrain the magnitude of the plane- 
tary field. These calculations assume that the accelerated 
electrons, which are the source for the radio signal, ap- 
pear during the period of time when the planetary and 
stellar fields are exactly opposite (in particular, during a 
geomagnetic event, when these fields can be opposite for 
several days). Our simulation shows that for the general 
case the situation of opposite fields is likely to be much 
shorter in the case of a realistic, more complex stellar 



G 



field, and that reconnection events last for an hour or 
two at the most. If radio emissions from such events can 
be detected at all, they will require long-term monitor- 
ing observations in order to distinguish such events from 
flares and other background radio signals. 

Despite its short duration, the reconnection event leads 
to a release of an amount of plasma from the planetary 
magnetosphere (the exact amount, of course, depends 
on the particular set of planetary boundary conditions). 
Based on the mass flux displayed in Figure[7j we define an 
average value of F = 1CP 11 g s^ 1 cmT 2 and a flux chan- 
nel cross section of A = ir(0AR+) 2 to determine the total 
planetary mass flux, M p = F ■ A = 1.5 • 1CT 14 Mj yr^ 1 . 
The average mass flux for the "Weak" case is about two 
times lower than the "Strong" case. If we assume that 
such reconnection events occur every orbit (2.2 days), the 
planetary magnetosphere loses 2.5 • 10~ 12 Mj in a single 
year (~ 2.5 • 10~ 15 M yr~ x ) — not a significant amount 
integrated over the lifetime of the planet. The recon- 
nection depends on the structure of the stellar magnetic 
field, which is not constant. Therefore, there are likely to 
be periods during which more reconnection events occur, 
and periods with fewer events. We further stress that this 
mass flux is of "cold" plasma escaping from the planetary 
magnetosphere and is not an attempt to model the actual 
atmospheric escape from the planet. Nevertheless, these 
reconnection events should include particle acceleration 
and precipitation, as well as radiation enhancement that 
might impact the atmospheric escape in general. In ad- 
dition, stellar coronal mass ejections could drive similar 
reconnection events on top of the reconnection between 
the ambient planetary and static stellar fields. Such ef- 
fects, however, are out of the scope of this paper and we 
defer further discussion to future work. 

It is tempting to associate the infalling planetary 
pl asma with the hot spot discovered in t he Ca II K band 
by |Shkolnik etaE] ( |2003| |2005a|b[ |20Q8[ >, but we caution 
that our simulations do not describe the precise mecha- 
nism by which energy is transferred from the planetary 
orbital motion to the stellar corona. The volumetric in- 
tegral of the sum of the magnetic and the kinetic energy 
in the simulation domain is 



5 • 10 4 s, so the total available power is: 



E,, 



FP 
8tt 



pu 

~2 



dV. 



(1) 



Throughout the simulation, the value of this integral for 
the "Strong" case ranges between (5.4 - 5.6)-10 3 ° ergs, 
with a change of the order of 1%. The change in the 
energy for the "Weak" case is about the same. The 
total volumetric energy for the MHD solution without 
the planet is 5.2 ■ 10 36 ergs so that the contribution of 
the planetary orbital motion and the planetary magnetic 
field is about 2 • 10 35 ergs. Since the MHD model is 
limited in its description of magnetic reconnection, it 
can only produce the change in magnetic field topol- 
ogy, which is the source for modulations in the total 
energy. In reality however, magnetic reconnection be- 
tween the stellar and planetary fields will involve dis- 
sipation of this energy, as well as non-thermal particle 
acceleration. Therefore, the 1% difference in the total 
energy, E avai i = 5 • 10 34 ergs, is the available energy for 
these mechanisms. The timescale, r, for the reconnection 
event in our simulation is about a quarter of an orbit, or 



P^E n 



i/r 



10 30 ergs s~ 



(2) 



Assuming an efficiency of 1% in converting this power, 
our simulation finds at least 10 28 ergs s _1 is available 
to accelerate particles along field lines connecting the 
planet and the star. These accelerated particles then 
hit the chromosphere and cause the observed hot-spots. 
This amount of energy is comparable with the total X- 



ray power of HD 189733 (Pillitteri et al. 20101, and it 



is reasonable to believe that it can produce the observed 
modulations in the Ca II band, which are also estimated 
to involve energies of the order of 10 2 8 — 10 28 ergs s~ 4 
( |Shkolnik et al.|2005aJjHall et al.|2007| ). Variations in the 
total mass flux at different radial distances in the sim- 
ulation domain (see Figure [5| reveal that the planetary 
magnetosphere contributes to the total mass flux dur- 
ing the reconnection event. The total mass flux through 
the Alfven surface based on the simulation without the 
planet is 8.3 • 10 -14 Mq yr~ Y , which is comparable to, 
but slightly smaller than, the values from the simulation 
that includes the planet. This differs from the result ob- 
tained in Paper II for a synchronously revolving planet: 
in that case, the magnetosphere served only as an obsta- 
cle that reduces the amount of stellar wind that crosses 
the Alfven surface. Here however, the planetary magne- 
tosphcric shape can be changed by the dynamics of the 
system, it can contribute to the total mass flux, and it 
disrupts the wind only for a short duration of time. It 
seems that for the particular set of parameters of the 
HD 189733 system, the effect of the planet on the stellar 
wind mass loss rate is not that important. 

5. SUMMARY AND CONCLUSIONS 

We have performed a realistic MHD simulation of SPI 
using the parameters of the HD 189733 planetary sys- 
tem. The simulation is based on the observed stellar 
surface magnetic field and includes the stellar wind, time- 
dependent planetary orbital motion and planetary out- 
flow. We choose to study two cases. One with strong 
planetary magnetic field and high base density and one 
with weak planetary field and low base density. Based 
on the simulation results, our main findings are: 

1. The use of a realistic stellar field results in a highly 
structured stellar corona characterized by sectors 
with different plasma properties. 

2. As the planet orbits the star, it interacts with 
the different sectors. Each sector crossing, as well 
as the transitions between the sectors, leads to 
changes in the SPI. 

3. The planet is followed by a long, comet-like tail 
of plasma that represents its magnetospheric tail 
wrapped around due to the rapid planetary orbital 
motion. This tail is almost perpendicular to the 
stellar wind direction and can lead to magneto- 
spheric dynamics which differ to those in the case 
in which the magnetospheric tail is aligned with the 
wind direction. 

4. A large reconnection event occurs when the plan- 
etary and stellar field are opposite in polarity, 
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where the reconnection site varies depends on the 
planetary field strength. This reconnection event 
causes a release of cold plasma from the plane- 
tary magnetosphere with a mass flux of M p — 
2.5 • 10 -12 Mj yr^ 1 . However, this escape of ma- 
terial can be defined as a magnetospheric escape 
rather than atmospheric escape. The latter cannot 
be define by the simulation presented here since it 
does not include a self-consistent atmospheric es- 
cape mechanism and is only controlled by a bound- 
ary value for the planetary density; 

5. The short duration of the reconnection event sug- 
gests that such SPI cannot be observed persistently, 
and in particular when such signatures are being 
searched for in radio bands. 

6. The simulation indicates that a sufficient amount 
of energy can be dissipated by the star-planet in- 
teraction to explain the hot-spots observed in the 
Ca II lines of stars hosting close-in planets. 

7. The total mass flux carried by the stellar wind is 



modulated by the planetary orbit and the partic- 
ular value at any time depends on the magnetic 
topology. 

In this work, we used the parameters of HD 189733, 
with planetary semi-major axis of a = 8.8i?*. It is rea- 
sonable to believe that with a smaller value of a, some 
of the SPI effects in the simulation would be enhanced. 
Further studies of SPI involving a physical atmospheric 
escape from the planet, as well as the impact of stellar 
coronal mass ejections on the planetary magnetosphere, 
would be highly motivated. 

We would like to thank Scott Wolk and Steve Saar 
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Space Weather Modeling Framework, developed by the 
Center for Space Environment Modeling, at the Univer- 
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NASA ESTO-CT, NSF KDI, and DoD MURI. 
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TABLE 1 

Stellar and Planetary Parameters of HD 189733 



Stellar Parameter Value 
R~i, 0.76 i? 
Mi, 0.82 M Q 
P* 11.95 d 

Planetary Parameter Value a A value of 2 has been adopted in this simulation) 

Rp 1-1 Rj a 
P orb 2.2 d 

B eq 2.0 G 

a 8.8 R* 




0.00 

Fig. 1. — A map of the radial stellar magnetic field reproduced from [Fares et al~] p010) in a polar projection (left), and in a longitude- 
latitude projection (right). 




Fig. 2. — Left: steady-state MHD solution driven by the surface field map without the planet. The stellar surface is shown as a 
sphere colored with contours of the radial magnetic field. Blue lines represent closed field lines and yellow lines represent open field lines, 
number density (middle) and B z (right) color contours are displayed on the equatorial plane along with magnetic filed lines shown as black 
streamlines. Selected phase angles of the planetary orbit are shown as well. 




Fig. 3. — Left: the grid structure adopted for the simulation including the planet displayed on the equatorial plane in black, with the 
star and the planet shown as red circles. Right: a zoom-in on the vicinity of the planet. 
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Fig. 4. — Color contours of the number density shown on the equatorial plane for selected phase angles for the "Strong" case. Axes are 
expressed in units of stellar radii. 
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Fig. 5. — Mass flux as a function of planetary phase angle integrated over spheres of r = 6i?* (blue), r = 10-R* (red), r = 23-R* (black). 
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Fig. 6. — Same display as Figure [4] but for the "Weak" case. Note that the minimum density value here is slightly lower than Figure [1] 
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Fig. 7. — Color contours of number density (left) and mass flux (right) displayed on the equatorial plane and on a sphere of r = 
with the planet shown as a blue sphere. White line separates positive and negative values of the mass flux. The plots are for the dynamic 
simulation for the time when a magnetic reconnection event has been identified, at phase angle of 60 degrees (top), and for the idealized 
simulation at a similar phase with an enhanced planetary surface density and dipolar stellar field (bottom). 
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Fig. 8. — Left four panels: the equatorial plane colored with number density contours, together with red and blue spheres representing the 
star and the planet, respectively. The panels are for phase angles of 270 (top-left), 60 (top-right), 100 (bottom- left), and 180 (bottom-right) 
degrees. Selected magnetic field lines are shown in yellow and blue. Two right panels: the equatorial plane colored with number density 
contours for phase angle of 60 degrees(top) and 100 degrees (bottom). The white arrows indicate the regions where reconnection is taking 
place. 
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Fig. 9. — Magnetospheric boundary at different phase angles defined by vanishing radial stellar wind. In each case the star lies to the 
right, out of the frame. The boundary is shown as a solid red line displayed on the star-planet meridional plain along with color contours 
of number density. Selected magnetic field lines are shown in white. Illustrated are phase angles of 270 (top-left panel), 180 (top-right 
panel), 60 (bottom-left panel), and 100 (bottom-right panel) degrees. 



